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ABSTRACT 



A geometrically similar family of three supercavitating 
hydrofoils was tested in the M.I.T. variable pressure water 
tunnel. The half-span foils were of elliptical planform; 
the ratio of foil half-span to tunnel height was 1/4 for the 
small foil, 1/2 for the medium foil, and 3/4 for the 
largest foil. The tunnel was of square cross section. Lift, 
drag, moment, tunnel speed, ambient pressure, and cavity 
pressure were measured for attack angles from 8 to 21 degrees 
and a variety of ambient pressure settings; cavity length 
measurements were obtained from photographs. Results were 
compared with theoretical results obtained by Leehey, and 
also with a more detailed numerical lifting surface procedure 
developed by Jiang and Leehey . For the small and medium foils, 
it was sufficient to correct only for the effect on downwash 
of the images of the trailing vortices. The large foil data, 
however, required further correction; upon application of the 
same corrections which were applied to the data for the two 
smaller foils, the force and moment data for the large foil 
plotted slightly lower than did the data for the two smaller 
foils, while the cavity length data for the large foil 
indicated cavity lengths significantly larger than for either 
of the theoretical predictions or the cavity length data for 
the two smaller foils. Through the application of existing 
two-dimensional corrections, the force data for the large 
foil were brought into close agreement with the force data 
for the two smaller foils, but no suitable correction factors 
exist for the cavity length data. 
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Title: Professor of Naval Architecture and 

Applied Mechanics 



-3- 



ACKNOWLEDGMENTS 



This research was carried out under the Naval Sea 
Systems Commcind, General Hydromechanics Research Program, 
Subproject SR 009 01 01, administered by the David W. Taylor 
Naval Ship Research and Development Center, Contract 
Number N00014-76-C-0358. The author is exceedingly grateful 
to Professor Patrick Leehey of the Department of Ocean 
Engineering at the Massachusetts Institute of Technology 
for his inspiration, suggestions, and guidance. Special 
thanks are also due to S. Dean Lewis and C. W. Jiang, also 
from M.I.T., who rendered valuable assistance in the 
performance of the experiments, and to Miss Cheryl Gibson 
for typing the various drafts and for enduring countless 
errors. Finally, the author is indebted to his wife, 
Gretchen, for her unfailing support and encouragement 
and for the many sacrifices endured by her during the 
preparation of this report. 



-4- 



TABLE OF CONTENTS 

Page 

TITLE PAGE 1 

ABSTRACT 2 

ACKNOWLEDGMENTS 3 

TABLE OF CONTENTS 4 

LIST OF FIGURES 5 

NOMENCLATURE 7 

INTRODUCTION 10 

EXPERIMENTAL APPARATUS AND PROCEDURE 12 

Apparatus 12 

Procedure 14 

RESULTS 17 

Lift 18 

Drag 19 

Moment 19 

Cavity Length 20 

Cavitation Number 21 

Influence of Static Pressure Controller 23 

DISCUSSION 26 

CONCLUSIONS 34 

REFERENCES 36 

APPENDIX I 73 

APPENDIX II 75 

APPENDIX III 77 

APPENDIX IV 81 



-5- 



LIST OF FIGURES 



Figure Page 

1 M.I.T. Test Facility 38 

2 Foil Dimensions 39 

3 Static Pressure Measurement 40 

4 Cj^/a^ vs. a^/a^, a = 8.0° 41 

5 a = 9.5° 42 

6 a = 11.0° 43 

7 a = 12.0° 44 

8 a = 14.0° 45 

9 a = 16.0° 46 

10 a = 18.0° 47 

11 a = 21.0° 48 

12 vs. a = 8.0° 49 

13 a = 9.5° 50 

14 a = 11.0° 51 

15 a = 12.0° 52 

16 a = 14.0° 53 

17 a = 16.0° 54 

18 a = 18.0° 55 

19 a = 21.0° 56 

20 

21 a = 9.5° 58 



22 a = 11.0° 



59 



- 6 - 



Figure Page 

23 ^ " 12.0° 60 

24 a = 14.0° 61 

25 a = 16.0° 62 

26 a = 18.0° 63 

27 a = 21.0° 64 

28 Wo, Small Foil 65 

29 Medium Foil 66 

30 Large Foil 67 

31 Comparison of Cavity Lengths 68 

32 vs. o„, Small Foil 69 

c V 

33 Medium Foil 70 

34 Large Foil 71 

35 Large Foil: Comparison of Cavity 

Pressure Measurement 
Techniques 



72 




m: 




-7- 





NOMENCLATURE 


a 


foil leading edge bevel length; see Fig. 2 


A 


2 

aspect ratio = 2s /S 


c 


foil mean chord, measured at centroid position 


c 


foil root chord 


6 


tunnel correction factor for influence of images 
of trailing vortex system 


s 


drag coefficient = -5^^ ~ ^ 

(l/2)pu2s ^uncor ^i 


Cd 

uncor 


uncorrected drag coefficient, as calculated from 
the raw experimental data = 

(l/ 2 )pu 2 s 




drag coefficient, corrected for blockage; see 
equation (1) 


1 


change in induced drag due to effect of images 
of trailing vortex system 




lift coefficient = (i/aipu^is 


a 

CM 

U 


two-dimensional lift coefficient; see equation 
(II-5) 




moment coefficient about the midchord 

. ^ moment 

point = = — 

(l/2)pu2sc 


H 


tunnel height = 51 cm 


^ (z) 


cavity length as a function of spanwise position 
measured from leading edge 


L 


cavity length at centroid position; measured 
from midchord 


Pb 


minimum pressure along tunnel wall in vicinity 
of cavity 
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Pc 



measured cavity pressure 



P 

Pu,D 




upstream tunnel static pressure at the spanwise 
position of the foil centroid 

vapor pressure of water 

tunnel static pressure controller setting 

static pressure signal to controller from 
upstream, downstream tunnel static pressure 
taps; see Fig. 3 

strength of discrete source located in the 
(i,j)-th element 



q(x,z) source strength at (x,z) 



s 



foil half-span 



As^j element length for (i,j)-th element 

S planform area of half-span model 



S tunnel test section cross-sectional area 

o 

t thickness of hydrofoil 



u,v perturbation velocities in the x,y-directions 

N 

U water velocity one meter upstream of dynamometer 

shaft axis 



V maximum water velocity (at point of minimum 

pressure, Pj^) in vicinity of cavity 

Xj^(z) chordwise coordinate of leading edge as a 

function of spanwise position 

distance of cavity pressure tap from midchord 

x,y,z foil coordinates; see Fig. 2 

distance of centroid of half-span foil model 
from midspan 

Z^ distance of cavity pressure tap from midspan 
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a 



a 



T 



Aa^ 



Y 

"ij 

Y (x,z) 
X 

^3D 






P 



o 



c 



o 



V 



o' 



a" 



geometric angle of attack with reference to 
tunnel side wall; when used in Appendices II 
and III, a indicates angle of attack in a free 
stream 



"true" angle of attack, as corrected for effects 
of images of trailing vortices = a + Aa^ 

change in angle of attack due to effect of 
images of trailing vortices 

arctan (2a/a^) 

strength of discrete vortex located in (i,j)-th 
element 



vortex strength at (x,z) 

ratio of foil frontal width to tunnel width; 
see equation (3) 

ratio of foil frontal area to tunnel cross- 
sectional area (without foil) ; see equation (4) 

dummy variables in the x, z-directions 

density of water 

cavitation number based on measured cavity. 

pressure = 

^ (1/2) PU2 



cavitation 
pressure = 



niiinber, based on vapor 

Poo-Pv 

(1/2)Pu2 



cavitation nxamber, corrected for blockage; see 
equations (2) and (5) 



cavitation number, based on 
the maximum water velocity = 



Pb' Pc' 
Ph-P c 
(1/2) pv2 
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INTRODUCTION 



Corrections for water tunnel wall effects during the 
testing of subcavitating hydrofoils are essentially the 
same corrections developed fifty years ago for use in wind 
tunnel testing of airfoils; references [l] and [2] are 
standard texts in this area. Only within the last twenty 
years have wall effects been considered in the case of 
supercavitating hydrofoils, with the majority of theoretical 
and experimental work concentrated on two-dimensional hydro- 
foils. The problem of blockage can be significant; the 
shape of the foil and cavity combination must be known in 
order to apply a blockage correction of the type utilized in 
subcavitating flow. Unfortunately, this proves to be an area 
of great difficulty in the various theories. Corrections to 
drag coefficient and cavitation number for the two-dimensional 
pure-drag case are given in references [3] through [7]; a 
nonlinear theoretical model of the two-dimensional lifting 
problem is formulated in reference [6]. In reference [8], 
Baker converts this model into a computer program which 
numerically calculates the wall effect on two-dimensional 
hydrofoils of arbitrary shape. It was hoped by Baker that 
these results would prove adequate for application to flows 
over high aspect ratio hydrofoils. Prior to the current 
research, a systematic series of experiments which exhibit 
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wall effects on a geometrically similar family of hydrofoils 
had not been carried out in either the two- or three- 
dimensional case. The purpose of the present research is 
to make available data on a geometrically similar family of 
three supercavitating hydrofoils of finite span, and to apply 
existing two-dimensional corrections in an attempt to bring 
the data into agreement with three-dimensional, unbounded 
flow theory. 
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EXPERIMENTAL APPARATUS AND PROCEDURE 



Apparatus 

The experiments were conducted in the M.I.T. recircula- 
ting, variable pressure water tunnel (Fig. 1) . The test 
section is 51 cm square in cross-section, and 147 cm long. 
Test section velocity was continuously variable between 0 
and 9.1 m/sec; a manometer indicated the difference in 
pressure across a contraction in the tunnel one meter up- 
stream of the dynamometer axis. This manometer had been 
previously calibrated using a pitot txibe in the test section. 
The assumption was made that the blockage caused by the foil 
and cavity combination did not affect the manometer 
appreciably. Test section static pressure was variable 
between 60 mm Hg absolute and atmospheric pressure; static 
pressure readings were taken using a mercury manometer 
connected to either or both of two taps on the centerline 
of a side wall of the test section. These taps were 
located 62 cm upstream and downstream of the hydrofoil 
dynamometer, which was located at the center of the test 
section. 

Foil dimensions are given in Fig. 2. The three half- 
span, aspect ratio five foils had sharp leading edges. 
Geometrically similar elliptical planforms were chosen even 
though foils of this type were more difficult to fabricate 
than, for example, foils of rectangular planform; elliptical 
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planform foils exhibit less tip loading effect than do 
rectangular planforin foils, thereby rendering them more 
amenable to analysis. Each foil was drilled with a hole 
terminating in a pressure tap on the cavity side of the 
foil. The medium-sized foil was the aspect ratio five foil 
utilized by Leehey and Stellinger [9], except that the 
cavity pressure tap was relocated to a position nearer the 
tip. The composition of the medixam foil was of type 304 
stainless steel, while that of the small and large foils 
was of type 614 alviminxim bronze. Each foil was welded to 
its own 10.16 cm diameter mounting plate; the mounting 
plates, in turn, were bolted to the dynamometer so that 
the small, medixim, and large foils protruded 1/4, 1/2, and 
3/4 of the way, respectively, downward from the top window 
into the test section. 

The dynamometer was capable of measuring three 
moment and three force components, but only the moment 
about midchord and the forces tangential to and perpendicular 
to the span and chord were measured. Zero geometrical 
angle of attack was set by aligning the flat (wetted) side 
of each foil with the test section side wall. Upon 
alignment, a sighting telescope, which was attached to the 
dynamometer, was aligned with a scale on the laboratory 
wall; this scale, graduated in increments of 0.05 degree, 
was utilized in setting the geometrical angle of attack. 
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Procedure 

Before test runs were begun on any one day, the tunnel 
was operated at medium velocity for twenty minutes at a 
pressure of 300 mm Hg absolute; this removed almost all of 
the air trapped in remote locations of the tunnel. Prior 
to each test run, the bubble chamber shown in Fig. 3 was 
utilized to collect and draw off any gases or air which may 
have been trapped in the static pressure sensing line; once 
cleared of gases, the bubble chamber was merely a reservoir 
of common pressure. On the other hand, small amounts of air 
were bled intermittently into the cavity pressure sensing 
line during testing in order to keep it free of moisture. 

Each foil was tested at a variety of attack angles 
from 8 to 21 degrees. The lower limit of 8 degrees was 
chosen so that the cavities would not terminate on the foil. 
For all test runs, the test section velocity was set at 
8.5 m/sec, whereupon static pressure was reduced until a 
cavity covered the entire planform. Readings were taken 
of lift, drag, moment, static pressure (with pressure taps 
U and D open; see Fig. 3) , cavity pressure, and velocity. 

A Polaroid photograph was taken of the foil and cavity from 
the wetted side of the foil, utilizing an exposure time 
of 1/40 second with bottom flood lighting. Pressure tap D 
was then closed, and the same readings taken; reopening 
pressure tap D and closing pressure tap U, the same series 
of readings was repeated once more. (Photographs were not 
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taken for these last two pressure tap combinations since 
hardly einy variation in cavity shape or length was noted 
when compared to the situation where both pressure taps 
were open . ) This completed the measurement of three data 
points at one controller setting. Both pressure taps 
were then opened, and the static pressure was reduced by 
20 to 30 mm Hg, whereupon the above procedure was repeated. 
The static pressure was reduced in this fashion until an 
excessive number of recirculated cavitation bubbles resulted 
in either extremely poor conditions for photography or 
rapidly fluctuating instrument readings. Prior to and 
immediately following a run of this type, the room and 
water temperatures were recorded, as were the tares on the 
instrumentation . 

A computer program was utilized in the data reduction. 
Wetted surface frictional drag (for each foil and its 
mounting plate) and dynamometer shaft twist were taken into 
account, and test section static pressure readings were 
modified for static head difference to give the static 
pressure at the vertical position of the foil centroid. 
(Gravity effects were otherwise neglected.) Cavitation 
number was computed using measured cavity pressure, and 
also by using vapor pressure at the ambient water temper- 
ature (computed using equation (2) on page 151 of 
reference [lO]) . Lift, drag, and moment coefficients 
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were also calculated. Standard wind tunnel corrections 
(as described in Appendix I) were applied to the data to 
account for the effect of the images of the trailing 
vortices; this resulted in an increase in drag coefficient 
and also in an increase in effective angle of attack. 
Utilizing the photographs, cavity lengths were measured 
from the midchord at the centroid position. 



RESULTS 



Experimental results were compared with predictions 
from Leehey's linearized theory [ll], which utilizes the 
method of matched asymptotic expansions; the theory is valid 
to first order in angle of attack, and to second order in 
the reciprocal of aspect ratio. Comparisons were also 
made with a more detailed linear lifting surface calculation 
developed by Jiang and Leehey [12], This numerical theory 
is expected to be more accurate at the lower values of 
Oc/oir^r where the cavity length is the same size as the foil 
span or longer; it should also be more accurate, in general, 
for prediction of moment coefficients. Synopses of the 
theories developed in references [ll] and [12] are given 
in Appendices II and III, respectively; the experimental 
data are given in Appendix IV, both in raw form and with 
standard wind tunnel corrections (Appendix I) applied. 

The wall boundary layer momentum thickness was 
approximately 1.9 mm, as inferred from measurements by 
Stellinger under similar conditions [9]. Since less than 
two percent of the wetted surface area of the small foil 
was within 1.9 mm of the upper wall, the upper wall 
boundary layer was assumed to have negligible effect on the 
force and moment coefficients of all three of the foils 



tested. 
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It should be noted that the only correction for wall 
effect which has been applied to the data in Figs. 4 
through 30 is the correction for the effects of the images 
of the trailing vortices (Appendix I) . Also worthy of note 
is the fact that the supercavitating condition was never 
achieved for the small foil at an angle of attack of 8 
degrees. This fact bears upon the lack of scaling evident 
in the force and moment data of Figs. 4, 12, and 20 for this 
foil. This matter will be considered again in the sub- 
sequent discussion section. 

Lift (Figs. 4 through 11 ) 

There was very good agreement, overall, between 
theoretical predictions and the experimental data. For 
small values of the similarity parameter (i.e., for 

long cavities) , there was much better agreement with the 
nvunerical theory than with the asymptotic theory; in this 
case, the foil and cavity combination no longer has a 
large aspect ratio. For the higher values of 
Leehey's theory appears to underpredict, especially at the 
higher angles of attack. This is in apparent contradiction 
with Fig. 4 of reference [9]; it should, however, be noted 
that the data of reference [9] were based on cavitation 
numbers calculated with vapor pressure rather than with 
measured cavity pressure. It can be seen from Figs. 32 
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through 35 of the present work that, when based on measured 
cavity pressure, the cavitation number will be smaller, so 
that the above contradiction is only apparent. The small 
and mediim foil data plotted almost on top of one another, 
whereas the large foil data plotted lower than did the small 
and mediim foil data, especially at attack angles greater 
than 11 degrees. 

Drag (Figs. 12 through 19 ) 

When the effects of streamwise foil or flow curvature 

are negligible, linear theory predicts so that 

2 2 

the drag data are plotted as Cj^/a^ vs. a^/a,p. plotted 

somewhat higher than the theoretical predictions except 
for very long or for very short cavities. As with the lift 
data, the small and medium foil data plotted together, with 
the large foil data plotting somewhat below these, the 
effect becoming more pronounced at the higher angles of 
attack . 

Moment (Figs. 20 through 27 ) 

The moment coefficient was taken about the midchord 
consistent with the right-hand rule. As pointed out by 
Leehey and Stellinger [9], Leehey's matched asymptotic 
expansion theory neglects lifting surface corrections 
(third order in the reciprocal of aspect ratio) for the 
moment coefficient. It is not surprising, therefore, that 
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the current experimental data agreed best with the more 
detailed lifting surface theory of Jiang and Leehey [12]. 

The load cell utilized for moment measurements on the small 
foil was of inadequate sensitivity, and consequently, much 
of the small foil data at the smaller attack angles may be 
less accurate than the data for the medium and large foils; 
comparisons must therefore be made between the large and the 
medivun foil data. Once again, the large foil data plotted 
below that of the medium foil, especially at the larger 
attack angles. 

Cavity Length (Figs. 28 through 31 ) 

As mentioned previously, cavity lengths were measured 
from midchord at the centroid position rather than from the 
leading edge, since this convention was utilized by Leehey 
Lll]. Cavity length (L) was nondimensionalized on the foil 
mean chord (c) . The cavity length data for the small and 
medium foils plotted on top of one another (Figs. 28 and 29); 
in both cases, the cavity lengths were slightly less than 
those predicted by theory, but the overall agreement with 
theory was good. In contrast, there was a marked deviation 
in the large foil cavity length (Fig. 30) and shape (Fig. 31), 
especially at the root, where local blockage was the worst. 
The tip vortices seen in the photographs had no evident 
effect upon the force or moment coefficients. Neither 
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Leehey's matched asymptotic expansion theory nor the numer- 
ical lifting surface theory developed by Jiang and Leehey 
accounts for the roll-up effects which produced these 
isolated trailing vortices in the wake. 

Representative photographs of the large, medium, and 
small foils at the same angle of attack and at approximately 
the same cavitation number are shown in Fig. 31. The cavity 
shapes shown were averaged over the exposure time of 1/40 
second; visual observations showed that the instantaneous 
cavity shape was highly unsteady. For the small and medium 
foils, the general outline of the cavity was elliptical, 
as predicted by theory; the cavity length was undoubtedly 
reduced toward the tip due to the hydrostatic pressure 
gradient existing in the tunnel. 

Cavitation Number (Figs. 32 through 35 ) 

Below a cavitation niomber of approximately 0.4, there 
was almost no difference between cavitation number calculated 
with measured cavity pressure or with vapor pressure. It 
should be noted that in this region there exists a large 
variation in cavity length from very short to very long. 

This indicates that for even short cavities, where ram 
effects of the re-entrant jet from the cavity end have 
heretofore been thought to have a pronounced effect on 
measured cavity pressure, the actual cavity pressure was 
approximately equal to the vapor pressure of the water. 
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Above a cavitation number of 0.4, the disparity between 
and increased with increasing cavitation number; the 
data in this region, however, were only for large angles of 
attack and small cavity lengths, indicating that there may 
indeed have been some effect caused by the impinging of the. 
re-entrant jet on the cavity pressure tap. 

After the series of basic experiments was completed, 
ram effects were further investigated on the large foil 
by Jiang and Leehey, this time utilizing a total head tube 
for cavity pressure measurement. The L-shaped total head 
tube protruded downwards into the cavity from the upper 
tunnel wall so that it was parallel to the foil wetted 
surface and pointed towards the leading edge and away from 
the impinging re-entrant jet. Cavity pressure readings 
were taken with both the standard foil surface pressure tap 
previously utilized and the total head tube, and vapor 
pressure was calculated as was done previously; results are 
shown in Fig. 35. For the larger attack angles and shorter 
cavities (higher cavitation nvimbers) , use of the total head 
tube significantly reduced re-entrant jet effects; the 
result was that cavitation nxambers calculated with total head 
tube data were as much as 5 to 8 percent higher than 
cavitation mnnbers calculated with data obtained from the 
previously used foil surface tap. Although a similar check 
was not made on the small and medium foils, it can be 
assumed that the same general trend would occur. It should 
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be noted, however, that data points in Figs. 32 and 33 
exhibit a more pronounced downward "hook" at higher 
cavitation nvunbers than do the large foil data (Figs. 34 
and 35), a fact which accounts for the noticeable "hooks" 
in the lift and moment data at the higher angles of 
attack (see, for example. Figs. 10, 11, 26, and 27) . 

If the total head t\abe had been utilized throughout the 
entire series of experiments , these deviations from theory 
at high attack angles and short cavities would probably 
have been less. 

Influence of Static Pressure Controller 

As stated earlier, readings were taken for the various 
combinations of the two static pressure taps; this was 
done in order to determine the effect of blockage on 
pressure and velocity at the two different tap positions. 

It was anticipated that this variation in static pressure 
tap combination would have no effect on the forces and 
moments experienced by the foil. Due to the location within 
the system of the static pressure controller (Fig. 3) 

(which was inadvertently overlooked) , the forces and moments 
varied significantly when the combination of static pressure 
taps was changed. The explanation for this is as follows. 
Before data were taken, both static pressure taps were opened 
and the controller setting reduced to p; the actual static 
pressure in the tunnel was decreased by the control system 
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until the combined signal from both taps was p. When the 
instrument readings settled out, data were taken. The 
downstream static pressure tap was then closed, so that the 
only signal to the controller was p^, which was greater than 
the controller setting, p. The controller therefore lowered 
the actual static pressure in the tunnel until py equaled p. 
Consequently, the effective cavitation number (as seen by the 
foil) decreased, causing the cavity to grow. This reduced 
the effective camber of the foil and cavity combination; 
this loss of effective camber resulted in a reduction of the 
force and moment coefficients from their previous values. 

On the other hand, the measured cavitation number increased 
slightly due to a small decrease in measured cavity pressure; 
the presence of the controller prevented the measured static 
pressure from changing significantly. The opposite effect 
was observed with the upstream tap open and the downstream 
tap closed (i.e., an increase in force and moment coefficients, 
an increase in effective cavitation number, and a slight 
decrease in measured cavitation number when compared to the 
case when both pressure taps were open) ; there were only a 
few of these data points taken. 

Although the desired pressure readings were not obtained 
in the current experiment (i.e., the actual tunnel static 
pressure was different for each of the tap combinations) , 
the results do show the significance of the blockage. When 
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comparing the results obtained with the different tap 
combinations for the small and medium foils, there was very 
little change in the plotted data points; the increased 
blockage in the case of the large foil resulted in signi- 
ficantly different force and moment data. For the sake 
of clarity, the data for the different tap combinations are 
shown only for the large foil at an angle of attack of 
14 degrees. Unless otherwise specified, further remarks 
concerning foil data will be for the data obtained with both 
static pressure taps open. 

The general trend of the data agrees with the experi- 
mental results obtained by Kermeen (see page 37 of reference 
[13]) ; he observed that the longer the cavity in super- 
cavitating flow, the smaller the force coefficients. The 
same conclusion follows from the theoretical results of 
Leehey, as may be seen from Fig. 4 of reference [9]. 

At a later time, the system indicated by the dashed 
lines in Fig. 3 will be utilized to obtain data on the large 
foil; the present system would be operated with only tap U 
open to provide an indication of static pressure at "infinity" 
to the controller. This configuration will allow pressure 
measurements to be made without affecting the input signal 
to the static pressure controller. 
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DISCUSSION 

For the small and mediiam foils, it appears to be 
necessary to correct only for the effects of the images of 
the trailing vortices in accordance with standard wind 
tunnel procedure. These corrections to drag coefficient 
and effective cuigle of attack brought the data into good 
agreement with theory; the lift, moment, and cavity length 
data also agreed well with theoretical predictions. For 
the most part, the plots for the small and medium foils 
were very close to each other. 

Upon application of the same corrections to the large 
foil data, good agreement with theory was seen for the force 
and moment coefficients; they were, however, generally 
lower than the data for the smaller foils, especially at 
higher attack angles. The plot of cavity length for the 
large foil, however, showed substantially longer cavities 
than predicted by either of the two theoretical models used 
for comparison. The most pronounced deviations in cavity 
length occurred for the long cavity (lower cr^/a,j,) data. 

When compared with the data obtained from the two 
smaller foils, the general trend of the large foil data is 
in agreement with Baker's predictions for wall effects on 
two-dimensional supercavitating hydrofoils [8]. Baker makes 
the point that the wall effect is negligible for tunnel 
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height-to-foil chord ratios (H/c) greater than 10 for 
cambered foils . Although he makes no such generalizations 
concerning supercavitating flat plates, (i.e., the case 
applicable to the current research) , it is seen from 
Figs. 4 through 27 that the effect of blockage on the force 
and moment coefficients is not exceedingly large. The 
following line of reasoning explains why the force and 
moment data for the large foil were somewhat lower than for 
the two smaller foils. For the three foils at the same 
cavitation nuinber and at the same angle of attack, the 
magnitude of the velocity on the cavity wall of all three 
foils is 0/2 (from linearized theory) ; consequently, all 
three foils have identical pressure coefficients on the 
cavity. The largest foil experiences significantly more 
blockage than does either of the smaller foils, so that the 
velocity on the wetted surface of the large foil is much 
higher than for either of the smaller foils; the pressure 
coefficient on the wetted surface of the largest foil is 
therefore the least of all three foils, resulting in smaller 
lift, drag, and moment coefficients than for either of the 
smaller foils. It is interesting to note that this overall 
effect is contrary to blockage effects experienced in 
subcavitating water tunnel or wind txinnel testing, where 
force and moment coefficients are increased. 
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Baker also questions the sharp increase in wall effect 
predicted by Wu, Whitney, and Lin [6] for thin symmetric 
wedges or for small attack angles in the lifting flow case. 
The results of the present work do not show a marked 
increase in wall effect on force and moment coefficients at 
the lower angles of attack, in agreement with Baker. 

Conversely, Baker's results indicate that the wall 
effect on the cavity plus wake can be tremendous. In the 
current research, this is, in fact, the most significant 
departure of the large foil data from the theoretical pre- 
dictions and from the data of the two smaller foils (Figs. 28 
through 30) . 

Several attempts were made to bring the large foil data 
into closer agreement with the small and medium foil data by 
using existing two-dimensional blockage corrections. Wu, 
Whitney, and Brennen [4] derive wall effect corrections 
for the two-dimensional pure-drag (wedge) case, employing 
the open-wake and Riabouchinsky models. Wu, Whitney, and 
Lin [6] consider wall effects in two-dimensional lifting 
cavity flows, a situation more applicable to the current 
research; Baker [8] utilizes a computer program to translate 
this into information which could be readily applied to 
correct experimental data. 
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The pure-drag, open-wake model employs corrections to 
both drag coefficient and cavitation number. The drag 
coefficient corrected for blockage, C^' , can be expressed as 



C ' = 



1+a' 



1+a. 



Cp + 0(X ) , 



( 1 ) 



where o', the cavitation number corrected for blockage, is 
given by 



1+0 



CpA + 0(X") , 



(2) 



where o^ is the measured cavitation number, is the drag 
coefficient as defined in the nomenclature, and X is the 
ratio of foil frontal width to tunnel width, which may be 
expressed as 

(mean chord) x sin (a) 

X= . (3) 

tunnel width 

For the large foil with X defined in this fashion, X = 0.30 

sin (a) . Application of the above corrections to the large 

foil data shifted the data points down and to the left on 

2 

the plots of Cp/a^ vs. a^/a^ (Figs. 12 through 19) but did 
not bring the data into any closer agreement with the small 
and medium foil data. When the cavity length data for the 
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large foil were replotted using the corrected cavitation 
numbers, the data points in Fig. 30 were moved downward; 
this downward movement was not sufficient, however, to bring 
the large foil cavity lengths into agreement with either 
theory or the small and meditim foil data. 

Actually, the straightforward application of such a 
two-dimensional correction to the case of a finite-span 
wing should overcorrect for blockage, making X much larger 
than it should be. Perhaps a better representation of X 
for the finite-span case would be 



foil frontal area 

tunnel cross-sectional area (without foil) 
(foil planform area) x sin (a) 
tunnel cross-sectional area (without foil) 



(4) 



Recalculating the corrections using place of X in 

equations (1) and (2) would make the corrections even smaller 
than previously calculated. Since the previous corrections 
proved to be inadequate, there is no reason to believe that 
the modified corrections should be any better. 

To correct measured cavitation niomber and measured 
drag coefficient utilizing the Riabouchinsky model, it is 
necessary to obtain pj^, the minimum pressure reading along 
the wall in the region of the cavity. Equation (1) is again 
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used to correct the drag coefficient, but now the corrected 
cavitation nvimber is calculated as 

a' = (2/3)a^ + (l/3)o' ' (5) 

where a'' is the cavitation number based on p^, pj^, and 
maximum velocity, V, and where is the cavitation number 
based on p^, p^, and upstream velocity, U. It had orginally 
been hoped that a reasonably close measurement of the minimum 
pressure could be obtained by using only the downstream tap 
(tap D in Fig. 3) ; as was shown earlier, the interaction of 
the static pressure controller precluded making this measure- 
ment accurately. In an attempt to apply the correction based 
on the Riabouchinsky model, the gross assumption was made 
that the maximum frontal area of the cavity was approximately 
the same as the frontal area of the foil; utilizing a 
continuity argument and the steady-flow Bernoulli equation, 
the pressure minimum was estimated, thus allowing application 
of the wall effect correction based on the Riabouchinsky 
model. Unfortianately , this procedure gave results which 
were comparedsle to those obtained using the open-wake model. 
The failure of both models in the current application 
should not be surprising since they were designed for two- 
dimensional, pure-drag cavity flows, and not for three- 
dimensional, lifting cavity flows. 
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Whereas the pure-drag wall effect corrections changed 
both drag coefficient and cavitation number, the results 
of Baker's computerization of the two-dimensional Wu, Whitney, 
and Lin [6] lifting case wall effect model are presented 
in terms of a change to the force coefficients (i.e., 
drag and lift coefficients in this case) , referenced to the 
uncorrected measured cavitation number, a^. Accordingly, 
Baker's corrections to the lift and drag coefficients are 
expressed as a percentage difference from the free stream 
condition. He presents data for only one ratio of tunnel 
height-to-arclength (i.e., chord length), H/c = 4. 

Calculating this ratio using the mean chord of the large foil 
gives H/c = 3.33; recalculating H/c as the ratio of tunnel 
cross-sectional area to foil area, H/c = 4.48 is obtained. 

It was concluded, therefore, that the H/c = 4 corrections 
given by Baker should be suitable for correction of the large 
foil data. Upon application of these additive corrections 
to both lift and drag data, it was found to bring the large 
foil data into close agreement with the small and medium 
foil data. 

Because an infinite wake model is utilized in 
reference [8], Baker's results give no correction for the 
extremely large cavity lengths recorded for the large foil. 

His results do show, however, that as blockage increases. 
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the nondimens ional cavity thickness increases. This was 
borne out nicely when the cavitation characteristics of the 
three foils utilized in the current research were compared 
at an angle of attack of 8 degrees. On the small foil, 
which produced the least blockage, supercavitation was never 
totally achieved, no matter how low the static pressure was 
brought. For the medium and large foils, which produced more 
blockage than did the small foil, cavity termination was 
achieved behind the foil over a range of static pressures. 
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CONCLUSIONS 



As a general rule, the agreement of the experimental 
data with unbounded flow theory was good for the small and 
medium foils, whose ratios of half-span to tunnel width 
were 1/4 and 1/2, respectively, and whose ratios of tunnel 
height to mean chord were approximately 10 and 5, respectively. 
It appeared that in order to bring the small foil and medium 
foil data into agreement with theory, it was necessary to 
apply only standard wind tunnel procedures in correcting 
for the effects of the images of the trailing vortices. 

Of the three geometrically similar aspect ratio five foils, 
the largest foil, which protruded 3/4 of the way into the 
tunnel, and which had a ratio of tunnel height to mean chord 
of approximately 3.33, had force and moment data which were 
close to the small and medium foil data, although slightly 
lower; the cavity length data for the large foil differed 
greatly from the small and medium foil data, showing much 
longer cavities for the same ratio of cavitation number to 
angle of attack. It is evident that blockage corrections 
are required in addition to the standard wind tunnel down- 
wash corrections. 

Application and modification of blockage corrections 
based on various two-dimensional pure-drag cavity flow models 
failed to bring the large foil data into agreement with the 
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data for the two smaller foils. It was found, however, that 
corrections based on a two-dimensional lifting cavity flow 
model brought the large foil force coefficient data into 
much closer agreement with the small and medium foil 
data. While these force coefficient corrections were found 
to give an "engineering order of accuracy, " no corrections 
have been found which give adequate corrections to the large 
foil cavity length data; it is clearly evident that further 
analytical work is necessary in this area. 

It is recommended that for future experiments, a static 
pressure measuring system similar to that depicted in Fig. 3 
be utilized. 
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APPENDIX I 

Standard Wind Tunnel Corrections for Effects of 
Images of Trailing Vortices 

The experimental data shown in Figs. 4 through 30 are 
corrected for the apparent upwash due to the interference 
of the tunnel walls, in accordance with standard wind 
tannel practice. The boundary condition of zero normal 
velocity at the tunnel walls is satisfied in the analytical 
treatment of the problem by introducing a doubly infinite 
system of image vortices; the interference due to the 
presence of the tunnel walls is then calculated to be the 
interference caused by the series of vortex images on the 
flow at the model. The net result is an additive correction 
to the induced drag and the angle of attack (to be perfectly 
correct, there is also a slight negative correction to the 
measured lift, but this is usually neglected) . Upon 
application of the correction factor, the flow about the 
model is transformed to the flow which would be seen if the 
foil were in an infinite stream. 

Since the hydrofoils tested in the square tunnel were 
only of half- span, a full- span model would require a duplex 
tunnel (that is, one whose ratio of height to width is 1/2). 
The ratio of model span to tunnel width for the small, medium, 
and large foils is 1/4, 1/2, and 3/4, respectively. Entering 
Fig. 6.34 of Pope and Harper [l], 6, the tunnel correction 
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factor for the influence of the images of the trailing 

vortex system was found to be 0.126 for the small, 0.104 for 

the medium, and Q.092 for the large foil. If the ratio of 

the foil area to txinnel cross-sectional area is S/S , then 

o 

the change in angle of attack is 

Aa. = 6^ Ct • (radians) (I-l) 

X O 

O 

The maximum value for Aa^ encountered in the experiments was 
1.47 degrees. When Aa^ is added to the measured geometrical 
angle of attack, a, the "true" angle of attack is obtained 

cxiji — u + Aot^ • (1—2) 



(In addition to the above correction, the data reduction 
program accounted for the dynamometer shaft twist, which 
never exceeded 0.06 degree.) The change in induced drag is 



AC 




C 



2 

L 



= Aa. 



(1-3) 



so that the total drag coefficient is 



= S 



+ AC 



uncor 



°i 



(1-4) 



where C_ is the uncorrected drag coefficient calculated 

^uncor 

from the raw experimental data. 
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APPENDIX II 

Leehey's Matched Asymptotic Solution for 
Supercavitating Hydrofoils of 
Large Aspect Ratio 



A synopsis of Leehey's theoretical results, as derived 
in references [9] and [ll], is given. 

Leehey's linearized steady-flow theory, which utilizes 
the method of matched asymptotic expansions, is valid to 
first order in angle of attack, and to second order in the 
reciprocal of the aspect ratio. Gravitational effects were 
neglected in the derivation. The analysis shows that the 
lift, drag, and moment coefficients, and non-dimensional 
cavity length (measured from midchord) are given by the 
following relationships: 

ira I 2sinY“l 



C. = 



siny ( 1+siny) 



1 - 



A (1+siny) 



(II-l) 



ira 






1 - 



siny (1+siny) 



2siny-l 1 
A (1+siny) j 



(II-2) 



-4a 



Cm - 



TT (1+siny) ^ 



1 - 



and 



l-2siny+2sin^y 
A siny (1+siny) 

8 



L/C = sec y - 1/2 - 



2 

a ' 



iraA 



'T ’ 

^2D 



(II-3) 



(II-4) 



where C- , the two-dimensional lift coefficient may be 
^2D 



expressed as 






(II-5) 



C 



L 



2D 



siny (1+siny) 



and where 

Y 



arctan 




and 



(1/2) pu2 ' 



(II-6) 



(II-7) 



Experiments performed by Leehey and Stellinger [9] showed 
that the above predictions for lift and drag coefficients 
and for cavity lengths were good; it was evident from the 
experimental results, however, that lifting surface corrections 
would be needed to improve the moment coefficient prediction. 



-77- 



APPENDIX III 

Numerical Lifting Surface Theory for Supercavitating 
Hydrofoils of Finite Span 

A brief description of the numerical lifting surface 
theory developed by Jiang and Leehey [12] for supercavitating 
hydrofoils of finite span will be given. 

The steady flow theory is linearized for application 
to thin wings at small attack angles. Discrete sources and 
vortices are utilized in the representation of the physical 
model, and the coupled integral equations are reduced to a 
set of simultaneous algebraic equations. The cavity closure 
condition is applied stripwise on the cavity in an iterative 
process to obtain the desired cavitation number over the 
supercavitated planform. 

This numerical theory should prove to be more accurate 
than the asymptotic theory of Leehey [ll]; an implicit 
assumption in the latter theory is that the cavity length is 
less than the span, so that for small o^/a (i.e., long 
cavities) , the numerical theory should indeed give better 
predictions. It should be noted, also, that for the 
a/a^ = 1.0 curves in Figs. 4, 5, 7, and 8 of reference [9], 
there is an ever-present difference between the experimental 
data and Leehey' s asymptotic theoretical predictions for 
lift and drag coefficients. For long and short cavities, 
however, there is an even stronger disagreement between data 
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and Leehey's predictions for moment coefficient (Figs. 9 
and 10 of reference [ 9 ]) . The numerical lifting surface 
theory was developed in an attempt to overcome these 
difficulties . 

The foil-cavity surface is collapsed to become a 
region in the (x-z) plane of the foil. The relationships 
for the jximp in the tangential and normal components of 
the perturbation velocity upon crossing the foil surface 
and cavity projection surface are: 

u(x,z,^0) ^ u(x,z, 0) = -y(x,z), (III-l) 

v(x,z,^0) - v(x,z, 0) = q(x,z), (III-2) 



where y and q are the vortex and source distributions, 
respectively. The boundary conditions are: 



v(x,z, 0) 



dy(x,z, 0) 



dx 



on the foil 
wetted surface 



u = a^/2' 



on the cavity 
boundary 



(III-4) 



q(C,C)dC = 0 
Jx^(z) 



cavity closure 
condition, (III-5) 



where Hz) is the cavity length, measured from the leading 
edge as a function of spanwise position, and Xj^(z) is the 
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chordwise coordinate of the leading edge as a function of 
spanwise position. 

The discrete vortex and source method is utilized in 
the formulation of a series of simultaneous equations. 

Upon stabdividing the projected surf ace into small elements, 
a discrete bound vortex, a trailing vortex, and a source 
are located within each of these elements. The bound 
vortex is situated on the quarter-chord line of each element, 
and induced velocities are calculated for all elements at 
their midspan, three-quarter chord positions. The 
concentrated source is taken to be a constant distribution 
across each element at its three-quarter chord position. 

The boundary condition of fixed cavity pressure is imposed 
on each element at its midspan, quarter-chord position. 
Utilizing these boundary conditions, the following equations 
are obtained; 
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(III-6) 



'kJl 



a 



1 Cl . . a 

+ _ y b hi - 0 

^ ^ h . °ijkJl a 

27T i,j 



a 



i = 1, . . . ,M (III-7) 
j = 1, . . . ,N 

= 0 for k > Ml 




j = 1, . . . ,N, (III-8) 



where and are the v and u velocities at the (k,il) 
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-th control point induced by the unit strength vortex and 
source at element (i,j). For the (i,j)-th element, is 

the element length, Yj^j the discrete vortex strength, and 
the discrete source strength. Ml and N are the number 
of elements on the foil along the chord and span, respectively. 
M is the number of elements along the chord in the cavity. 
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APPENDIX IV 

Data Reduction Program and Data Printout 

Beginning on p. 82, a listing is given of the FORTRAN IV 
computer program utilized in the reduction of the experimental 
data. 

The program output is also given, beginning on p. 92. 

For each foil at each angle of attack, a complete printout 
of the experimental data is given in raw form, with corrections 
for instrument tares, in reduced fomn, in nondimensional form 
without any corrections for wall effect, and finally in 
nondimensional form with corrections for the effects of the 
images of the trailing vortex system. 
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